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Abstract. The latter author, together with collaborators, pro- 
posed a numerical scheme to calculate the price of barrier options in 
[2l IH [5] . The scheme is based on a symmetrization of diffusion pro- 
cess. The present paper aims to give a mathematical credit to the 
use of the numerical scheme for Heston or SABR type stochastic 
volatility models. This will be done by showing a fairly general re- 
sult on the symmetrization (in multi-dimension/multi- reflections). 
Further applications (to time-inhomogeneous diffusions/ to time 
dependent boundaries/to curved boundaries) are also discussed. 

1. Introduction 

Recently, the latter author, together with Yuta Ishigaki, Takuya 
Kawagoe and Toshiki Okumura, introduced a new numerical scheme 
for calculating the price of barrier options in a series of papers [3j HI 
5J. The scheme is based on what they call "put-call symmetry", a 
notion introduced by Peter Carr and Roger Lee jT] in relation to (a 
generalization of) static hedging of barrier options. 

The put-call symmetry, PCS for short, of a diffusion process X at a 
point K e R, means the equivalence in law between K — X t (put) and 
X t — K (call) for arbitrary t > tk, where tk is the first hitting time at 
K. This is much weaker than the reflection principle which has been 
widely-recognized as a fundamental requirement for the static hedging 
of barrier options. While the put-call symmetry is still something one 
cannot expect for a given diffusion without luck, the latter author and 
her collaborators have noticed that, for any given diffusion 

• one can construct (easily) another diffusion that is identical 
with the diffusion up to the first hitting time and satisfies the 
put-call symmetry, 

• by which the price of barrier- type options written on the original 
diffusion is expressed by a combination of the prices of plain 
options written on the constructed one. 

• The fact in turn implies that the "symmetrization" offers a new 

numerical scheme for calculating the price of barrier options; it 
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transforms path-dependent expectations to path-independent 
ones. 

The numerical experiments they performed show that the scheme is 
quite plausible. They also claim in [U E] that it can be applicable 
to stochastic volatility models where the stock and its volatility are 
described by a two-dimensional diffusion. 

In this paper, we give a mathematical backgrounds for the scheme 
by establishing some symmetry results in a more general setting under 
the action of a reflection group. This in turn leads to further possible 
applications of the scheme. We will show that it can be used for the 
cases with time-inhomogeneous diffusions, time-dependent boundaries, 
as well as curved boundaries. 

This paper is organized as follows. We start with a detailed dis- 
cussion of the put-call symmetry (section |2J). After recalling the one- 
dimensional cases (subsection !2.ip . we first introduce a multi-dimensional 
generalization (subsection I2.2p and then extends it to multi-reflections 
(subsection 12. 3p . The key assumptions are found in the statement of 
Lemma 12.21 Section [3] is devoted to discussions on the 'symmetriza- 
tion" . Again staring from a review of one- dimensional cases (subsection 
I3.ip . we give the main result (Theorem I3.4[) in full generality in subsec- 
tion 13.21 Applications of the main theorem are presented in section HI 
We first give a credit to the use of the scheme for stochastic volatility 
models including Heston's and SABR type (subsection 14. ip . The trick 
for the generic stochastic volatility models is extended to get a general- 
ization of Theorem 13.41 (Theorem H~3"j) . as a corollary to which we show 
that the scheme is applicable even to time-inhomogeneous diffusions 
(Corollary 14.4(1 . This observation further enables the applications to 
curved boundary cases (subsection 14.31) and time-dependent boundary 
cases (subsection 14.4(1 . 

2. The put-call symmetry 

2.1. One dimensional case revisited. Let X be a one dimensional 
diffusion satisfying the put-call symmetry at K e R; 

(2.1) X t l {TK <t } = (2K - X t )l {TK < t} 

for any t > 0, where r K = inf{t > : X t = K}. The PCS (123]) 
alone, without the reflection principle, suffices to have a static hedging 
formula for barrier options. For completeness, we give a brief proof on 
it. Put 

X t = X t l{ TK>t j + (2K — X t )l{ TA -< t }. 
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Then by the PCS (JSHJ, we have 

(2.2) X t i X t; 

for arbitrary £ > 0. The equivalence in law (12. 2p in turn implies 

(2.3) P(X t eA:r K <t) = P(X t G A : r K < t) 

for any Borel set A. Suppose that Xq > K and A G {x > K}. Then 
the right-hand-side is equal to P{2K — X t G A) since {tk < £} is 
included in {2K — X t G A}. Thus it holds that 

(2.4) P(X t eA:r K >t) = P(X t G A) - P(X t G 2K — A) 

for any Borel set A with A G {x > K} (see pQ, [3]). In other words, 

E[f(X t -K)l {Xt>K} l {TK>t} ] 
1 ' } = E[f{X t - K)l {Xt>K} ) - E[f(K - X t )l {Xt<K} ] 

for any bounded Borel function / and t > 0, which can be understood 
as a static hedging formula. 

2.2. A multi-dimensional generalization of PCS. To generalize 
the argument in the previous subsection, we understand x i— > 2k — x 
as a reflection. In R d , a reflection is associated with a hyperplane. A 
hyperplane is given by Pl a ^ '■= {x G R d : (a,x) = k}, where k G R 
and 0) G R d . The reflection s at k '■ R d — > R d with respect to H a ^ 
is given by 

2a 

s a ,k(x) = x - ((x, a) - k)—-=. 

\ a \ 

Notice that s 2 a k = Id R d. A natural extension of the previous PCS could 
be as follows. 

Definition 2.1. Let X be a diffusion process in R d . We say X has 
the put-call symmetry with respect to the hyperplane H a j, if 

(2.6) X tl{T a>k <t} = SaAXt) 1 ^,^} 
for any t > 0, where r a ^ = inf{£ > : X t = Har- 
iri a totally similar way as above, we can obtain a static hedging 

formula corresponding to f)2.5p . In fact, we have, for t > 0, 

(2.7) E[f(X t )l {Ta k>t} ] = E[f(X t )} - E[f(s a , k (X t ))]. 

for any bounded measurable / whose support is included in a half space 
{x G R d | ±((a, x) — k) > 0}, provided that X has the PCS with respect 
to a hyperplane H a k . 
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2.3. The PCS with respect to a reflection group. As the reflec- 
tion principle is generalized to multiple reflections which form a group 
(see [6] and references therein), so is the put-call symmetry. In this 
section, we discuss the generalization in detail. 

It may be natural that we consider the exit time out of an intersection 
of hyperplanes. More precisely, denoting 

H+ k := {xeR d \(a,x)~k> 0} 

and 

:= p| if+ fc , 

(a,fc)e* 

we set, for a given diffusion X, 

r s<1> := mf{t > : X t £$} = min r a k , 

(a,fe)gf> 

and consider the problem of representing the expectation of /(JQ)l{ Ts >t } 
by those of f(g(X t )), where g runs through a set G, which will turn 
out to be the group generated by the reflections s\, A G $. 

Looking at the discussion in section I2.1[ we notice that the key was 
the equation (12. 3p . which is not anymore directly applicable to the 
multi-reflection case. However, we have the following generalization. 

Lemma 2.2. Let G be the group generated by the reflections {s Qj fc : 
(a, k) G $} and X be a diffusion process in H d satisfying PCS with 
respect to H ak for all (a,k) G $. Assume that (i) g£ D g'H = 
whenever g ^ g' G G, and (ii) there is a group homomorphism rj : G — > 
C (character) such that rj(s a ^) = —1 for each (a, k) G $. Then, for a 
Borel subset A ofH, x G £, and t > 0, we have 

(2.8) itiWXt G gA : r s# <t) = 0. 

geG 

Proof. We first note that the left-hand-side of (12.81) is absolutely con- 
vergent since the sets gA, g G G are disjoint. Therefore we can change 
the order as 

J2v(9)Px(X t egA:r^<t) 

(2-9) ^ _ _ 

= 12 z2 ^)Px{X t G gA : r E$ = r a>fc < t). 
(a,fc)e* g&G 

By the assumption on the put-call symmetry, we have 
(2.10) 

P x {X t G gA : r s$ = r a<k < t) = P x {s a>k (X t ) G gA : r s<6 = r Qjfc < t) 

= P x {X t G s aik gA : r s$ = r a ,* < £). 



On the other hand, by the assumption on rj we have 
} j rj(g)P x (X t G s a<k gA : r s<s = r a>fc < t) 

geG 

= v(s a ,kg)Px(X t G s a , k gA : r s$ = t QjA . < t), 

geG 

which is equal to 

- v(g)Px(X t G gA : r s$ = r a , fc < t), 
g eG 

thanks to the group structure. This observation together with the 
equation (I2.10p shows that 

J2v(9)Px(X t G gA : r E4> = r Qjfc < t) = 0, 

which proves the assertion with (12.91) . □ 

Thanks to the lemma, we have a generalization of the static hedging 
formula. 

Theorem 2.3. We keep the notations and the assumptions of the 
Lemma lKB . For a bounded Borel function f with its support in we 
have 

(2.11) E[f(X t )l { ^ >t} ) = X^WGT 1 ^))]. 

geG 

Proof. It suffices to show that 

P x (X t eA:r^>t) = J2v(g)P x (X t G gA), 

geG 

but this is equivalent to (12.81) since 

P x {X t G gA : r Stt > t) = 
unless g 7^ 1, by the assumption (ii) in Lemma [2. 2[ □ 

Remark 2.4. The one reflection case in section [2721 automatically sat- 
isfies the assumptions in Lemma I2.2[ and apparently (12. lip includes 

(fm 

Remark 2.5. The case with $ = {(a, fc), (—a, —A; — 1)} corresponds 
to the double boundary reflections, where, by choosing conventionally 
I a I = 1 and rj to be the (mod 2) "length" of the reflections (see e.g. [2]), 

P x (X t G A;t^ >t) = J2i P *( X t e A+2ka)-P x (X t G -A+(2k+l)a)}. 

kez 
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Remark 2.6. More generally, in the cases of multi-reflections, there 
do exist the situations where the assumptions are fulfilled. There are 
two important classes. One is that of finite reflection groups and the 
other, of affine reflection groups. The domain £ is a cone in the former 
class and a simplex in the latter. In fact, if {a} forms a so-called simple 
system (or fundamental system) of a root system, and (i) if k is fixed, 
then the group becomes finite, and (ii) if additionally a is taken from 
a root system properly and k' is taken to be —k — 1, then the group 
becomes (isomorphic) to the semi-direct product of the finite reflection 
group and the translation group Z, which is called an affine reflection 
group. In both cases we can take j] to be its determinant when the finite 
reflection group is embedded into the orthogonal group and therefore 
into the general linear group. For details, see e.g. [2]. 

3. The symmetrization 
Let X be a solution to 

(3.1) dX t = a{X t ) dW t + fi(X t ) dt, 

where a : R d — > H d x R d and /i : R d — > R d are piecewise continuous 
functions with at most linear growth, and W denotes a (/-dimensional 
standard Wiener process. We further impose some ellipticity condtions 
on era* to ensure that a unique (weak) solution to (13.11) exists, and that 
Euler-Maruyama approximation of the solution of (13.11) works; 

3.1. One dimensional case revisited. In the case of d = 1 (with 
a = 1), Carr and Lee [TJ showed that 

(3.2) a(x) = ±a(2k-x), ^(x) = ~n(2k - x) 

is a sufficient condition under which the solution X of (13. ip has the 
PCS. As is pointed our in Introduction, even though this is much 
weaker than the reflection principle which can be now rephrased as 
"PCS at any k", yet it is not practical to assume a price process to 
satisfy PCS. 

In [3], however, the following observation is made; 

Theorem 3.1 ([3J). Let X be a solution to \3. 1\) without \3. ty) and X 
be one to 

dX t = a(X t )dW t + fi(X t )dt, 

where 

a(x) = a(x)l {x>k} ± a(2k - x)l {x < k} , 

and 

fl(x) = n(x)l {x>k} - fi(2k - x)l {x < k} . 
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We assume Xq = Xq > K . Then we have 
E[f(X t -K)l {Xt>K} l {TK>t} ] 

= E[f(Xt - K)l {jCt>K} ] - E[f(K - X t )l {jtt<K} ] 
for any bounded Borel function f and t > 0. 

They claimed that the formula (13.31) gives a new insight to financial 
engineering of barrier options; it says that the price of barrier option 
is expressed by those of plain options. Numerical analysis for the for- 
mer is difficult, while the latter is much easier. In fact, the numerical 
experiments they performed (part of which is appeared in ([3])) show 
that their scheme is quiet effective. 

Remark 3.2. The PCS method can be seen as a diffusion equation 
counterpart (or simply a generalization) of so-called method of image 
charges in electrostatics. The authors thank Prof. Sergey Nadtochiy 
for suggesting to us this "synchronicity" . 

3.2. Symmetrization with respect to a reflection group. They 
also applied their scheme to stochastic volatility models like Heston or 
SABR type, where above Theorem 13.11 cannot be applied directly since 
they are basically two dimensional models. With a view to endow- 
ing a certificate, we give a general result on mult i- dimensional multi- 
reflection cases in this section. The certificate will be provided in the 
next section corollary to the result. 
We start with a lemma. 

Lemma 3.3. Let A be an affine transformation in R d such that Ax = 
Aqx + a for x G H d with Aq G GL(d, R) and a G R d . Suppose that 

(3.4) a(Ax) = A a(x)U x , fi(Ax) = A)/i(x) 

forx G R d , with some piecewise continuous x i->- U x G 0(d). Then AX t 
starting from Ax is identically distributed as X t starting from x G R d 
as a stochastic process provided that Ax = x, where X is the unique 
weak solution to A3. 

Proof. Put Y = AX. Then, 

dY t = A dX t = A a(X t )dW t + A of x(X t )dt, 
which equals to 

a{AX t ) U- l dW t + fi(AX t ) dt, 

by the assumptions (13 ,4p . where we note that U~ x W t is another Wiener 
process. Namely, Y is a weak solution to 

dY t = a(Y t )dW t + »(Y t )dt. 
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By the uniqueness of the solution, we have the assertion. □ 

As in the one dimensional case, we symmetrize both a and \i to get 
another diffusion with the PCS for which an extended static hedging 
relation still holds. 

Theorem 3.4. Suppose that we are given a family of hyperplanes in- 
dexed by $ which satisfies the assumptions of Lemma \2.2i Put 
(3.5) 

geG geG 

where T g is an orthogonal matrix corresponding to the reflection part 
of g, and x U x G 0(d) is a piecewise continuous map. Let X be a 
solution to 

dX t = a(X t ) dW t + j2(X t ) dt. 
Then for a bounded Borel function f with its support in we have 

(3.6) E[f(x t )i { ^ >t} ] = Y J v(g)E[fXg-\x t ))}. 

geG 

Here we used the notations set in section Wlh 

Proof. It suffices to show that a and fi satisfy f)3.4p . For h G G, we can 
find ho G R d and T/, G 0(d) such that hx = ThX + h^. We then have 

cr(hx) = ^ThT^TgCr^hx^^h-igx^Uhx 
gee 

= T h y^T h -i g a((h~ l gy 1 x)l {xeh -ig^ } U hx 

geG 

= T h a(x)U hx . 

Here we have used TfT g = Tj g for f,g £ G, which can be verified by 

fg(x) = f(T g x + g ) = T f TgX + T f g + f , 

where fx = TfX + fo and gx = T g x + go. 

Similarly we have fi(hx) = Thfi(x). □ 

Remark 3.5. We can use the numerical scheme based on the equation 
(13.61) for general $ since, by a linear transformation, a generic $ (and 
given X) can be transformed into the system satisfying the assumptions 
in Lemma 12.21 
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4. Applications 



4.1. Stochastic volatility models. As we have stated in the previous 
section, we give a certificate to the use of the PCS symmetrization 
method to such stochastic volatility models as Heston's and SABR 
type. 

A generic stochastic volatility model is given as follows: 

dX t = an(X t ,V t )dW t + fi 1 (X t ,V t )dt 

{ ' 1 ' dV t = a 2l (V t )dW t + a 22 {V t )dB t + fj, 2 (V t ) dt, 

where W and B are mutually independent (1-dim) Wiener processes, 
fi(x,v) = ()Ui(x, v), [i 2 {v )) is a continuous function on R 2 , and 

a(xv)= (°ntyv) 

is a continuous function on R 2 . In most cases, an(x,v) = xu(v) for 
some v and \h\{x,v) = rx from financial reasoning, and x- and v- axes 
are set to be natural boundary. In our scheme, however, these features 
are irrelevant but that V do not depend on S is important. 

Corollary 4.1. Let K > and put 

. \a n (x,v) x > K 

I -a n (2K - x,v) x<K 

/J,i(x, v) x > K 

-[ii(2K — x, v ) x<K" 

and let X be the unique (weak) solution to 

dX t = a n (X u V t )dW t + MXt, V t ) dt, 

where V is the solution to ^4-l\) - Then, it holds for any bounded Borel 
function f and t > that 

E[f(X t -K)l {Xt>K} l {TK>t} ] 

= E[f(X t - K)l {jCt>K} ] - E[f(K - X t )l { x t<K} ], 

where X is the solution to (f^.-Zp with X > K and tk is the first hitting 
time of X to K. 

Proof. Apply Theorem E31 to $ = {((1, 0), K)}, with U x = I. □ 

We can also obtain a static hedging formula for a double barrier 
option. Let K' be a positive constant. We consider static hedging of 
an option knocked out if X of (14. ip hit either the boundary x = K or 
x = K + K'. 



Corollary 4.2. Set 

a u (x,v) 

ja n (x-2nK',v) K + 2nK' < x < k + {2n + l)K', 

" \-a n (2K -x + 2nK',v) K +(2n-l)K' <x< K + 2nK', 

n G Z, 

p,i(x,v) 

//i(x - 2nK',v) K + 2nK' <x< k + (2n + l)K', 

-Hi(2K - x + 2nK', v) K + (2n - 1)K' < x < K + 2nK', 

n G Z, 

and let X be the unique (weak) solution to 

dX t = a n (X t , V t )dW t + fn(X t , V t ) dt, 

where V is the solution to h4-l\ ). Then, it holds for any bounded Borel 
function f and t > that 

E[f(X t - K)l { x t e(K,K+K')}l{r {KiK+Kl) >t}} 

= 5>[/(l t — K — 2nK')l {jtteiK>K+K , )} ) 

neZ 

- J2 E[f{K ~X t + (2n - 1)K')1^ MK>K+KI)} ], 

neZ 

where X is the solution to fl^.ip X G (K, K + K') and T(k,k+k') 
is the first exit time of X out of (K, K + K'). 

Proof. Apply Theorem [331 to $ = {((1, 0), K), ((-1, 0), K+K')}, with 
U x = 1. □ 

4.2. Time-inhomogeneous extension. The trick of Corollary 14.11 
further leads to a time-inhomogeneous extension of Theorem 13.41 Even 
more generally, we have the following 



Theorem 4.3. We reset X and V in (f^.ip as follows: X is R 1 valued 



and V is R d2 -valued, with an : R dl x FT 2 ->■ R dl g> R d \ a 2i : R d2 -> 
R d 2 R d % j = i ? 2, //! : R dl ® R d2 ->■ R d2 ; and ^ 2 : R d2 -> R d2 , all 
piecewise continuous, and W , d\- dimensional, and B, d%- dimensional, 
are mutually independent Wiener processes. For the SDE for V , we 
simply assume that = 1,2 and 1x2 are given so that there is a unique 
strong solution. Suppose that $ C R dl ® R satisfies the assumptions 
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of Lemma \2. 6 A Let G be the reflection group associated with $. Then 
we still have A3. 6]R 



Proof. Set $ = {((a, 0),k) G R dl+d2 : (a, k) G $ }, and apply Theorem 
1531 □ 

Corollary 4.4. Theorem \3.J\ is generalized to time-inhomogeneous cases 
o and fx can be time- dependent. 

Proof. Take d'i — \ and o~?\ = 022 = (and fi2 = 1) in Theorem 14. 31 □ 

4.3. Curved boundary. Let D C R d be a domain homeomorphic to 
a satisfying the assumptions in Lemma 12. 21 The scheme can be 
applied to the numerical calculation of such an expectation as 

E [f( X t)l{x t eD}l{T D >t}], 

where 

t d := mi{t > : X t £ D}. 

In fact, at least heuristically, F : D — > being a homeomorphism, 
Y t = F(X t ) can be realized as a solution to a time-inhomogeneous SDE 
in and we can apply Corollary 14.41 Note that in the scheme the 
information of the outside of of the coefficients of SDE is totally 
irrelevant. 

4.4. Time-dependent boundary. The scheme is also applicable to 
the cases where knock-out boundaries are time-dependent. 

Let an : [0, 00) — y R d , i — 1, • • • , d be smooth maps such that the 
matrix A(t) = [ai(t), ■ ■ ■ ,otd(t)] is invertible for all t and is at most 
linear growth in t. For each t > 0, put $(t) = {(aij(i), and let 
S(i) be the intersection of the half spaces associated with $(£). We 
will briefly discuss how it can be transformed to the previous problem, 
though heuristically. 

Let C : [0, 00) — > H d (g) R d be a solution (unique up to the initial 
point) to the following matrix valued linear differential equation: 

(4.2) C'(t) = C t A{t){A{t)Y l . 

Then C(t)ai(t) = (7(0)^(0) =: on for z = 1, 2, • • ■ , d since (g^D is 
equivalent to 

(C(t)A(t))> = 0. 

Moreover, since 

(a t (t),x) = (C^a^^^ty'x) = (a^C^ty'x), 

1 Precisely speaking, the group G should be embedded in Afhne group in H dl+d2 
2 Without loss of generality, we may assume that ki, i = 1, 2, • • • , d are constants. 
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we see that x t— > C*(t)~ 1 x transforms £(£) to a time-independent do- 
main £$» with $* = ki),i — 1, 2, • • • , d}. Now the problem turns 
into the pricing of options written on C*(t)~ 1 X t knocked out at the 
boundary of £$* . As we have commented in Remark 13. 5[ we can as- 
sume that $* satisfies the assumptions in Lemma l2~2l Combining them 
all, we see that Corollary 14.41 is now applicable. 
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